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Abstract 

We extend to incremental thermoelectroelasticity with biasing fields 
certain classical theorems, that have been stated and proved in linear 
thermopiezoelectricity referred to a natural configuration. A unique- 
ness theorem for the solutions to the initial boundary value problem, 
the generalized Hamilton principle and a theorem of reciprocity of 
work are deduced for incremental fields superposed on finite biasing 
fields in a thermoelectroelastic body. 

1 Introduction 

In the last decades, with the increasing wide use in sensing and actuation, 
materials exhibiting couplings between elastic, electric, magnetic and thermal 
fields have attrached much attention. 

In order to give certainty to experimental results and applications, the 
interest of many researchers turned to the mathematical fitting of these top- 
ics. 

Many applications have their mathematical formulation within a linear 
framework, and the theoretical study began from this context. 

Foundamental is Nowacki's paper [1], where a uniqueness theorem for 
the solutions of the initial boundary value problems is proved in linear ther- 
mopiezoelectricity referred to a natural state, i.e., without biasing (or initial) 
fields. Hence Nowacki j2] also deduced a generalized Hamilton principle and 
a theorem of reciprocity of work. 
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Li [3] generalized the uniqueness and reciprocity theorems for linear thermo- 
electro-magneto-elasticity referred to a natural state. 

Aouadi [I] establishes a reciprocal theorem for a linear theory in which 
the heat flux is considered as a constitutive independent variable, a rate-type 
evolution equation for it is added to the system of constitutive equations, and 
the entropy inequality is stated in the form proposed by Miiller [5]. 

Iesan [U] uses the Green-Naghdi theory of thermomechanics of continua 
to derive a linear theory of thermoelasticity with internal structure where in 
particular a uniqueness result holds. 

Related works on thermoelasticity and thermoelectromagnetism can be 
found in [ZJ to mj . 

The classical linear theory of thermopiezoelectricity assumes infinitesimal 
deviations of the field variables from a reference state, where there are no 
initial mechanical and electric fields. In order to describe the response of 
thermoelectroelastic materials in the presence of initial fields one needs the 
theory for infinitesimal fields superposed on initial fields, and this can only 
be derived from the fully nonlinear theory of thermoelectroelasticity. The 
equations of nonlinear thermoelectroelasticity were given in Tiersten |12j . 
Yang [13J then derived from [T2] the equations for infinitesimal incremental 
fields superposed on finite biasing fields in a thermoelectroelastic body with 
no assumption on the biasing fields. 

Here we extend the aforementioned three Nowacki's theorems pQ, [2] to 
incremental thermoelectroelasticity with initial fields. 

We explicitly refer to the incremental theory [13J, hence we rewrite from 
this paper, with the same notations, some formulae and results on constitu- 
tive equations of incremental thermoelectroelasticity. 

Of course, the theorems proved here just reduce to the ones in Nowacki's 
[2] by neglecting the initial fields. 

In the uniqueness theorem in Section H] we assume that in the initial 
state entropy does not depend on time and temperature is uniform. For the 
theorem of reciprocity of work in Section [6] we assume that in the initial state 
both entropy and temperature fields do not depend on time. 

2 Equations of Nonlinear Thermoelectroelas- 
ticity 

2.1 Balance laws and constitutive equations 

Consider a thermoelectroelastic body that, in the reference configuration, 
occupies a region V with boundary surface S. The motion of the body is 
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described by 

Hi = Vi(X L , t) , 

where denotes the present coordinates and X L the reference coordinates 
of material points with respect to the same Cartesian coordinate system. 

Let Kij, p , fj, A/,, pe, 0, 77, Ql and 7 respectively denote the first 
Piola-Kirchoff stress tensor, the mass density in the reference configuration, 
the body force per unit mass, the reference electric displacement vector, the 
free charge density per unit undeformed volume, the absolute temperature, 
the entropy per unit mass, the reference heat flux vector, and the body 
heat source per unit mass. Then we have the following equations of motion, 
electrostatics, and heat conduction written in material form with respect to 
the reference configuration: 

K L i,L + Pofi = PoVi, (1) 

Al,l = Pe , (2) 
PoOrj = -Ql,l + Pol , (3) 

The above equations are adjoined by constitutive relations defined by the 
specification of the free energy ip and heat flux Ql. 

V = i/>(E MN , W M , 9), Q L = Ql(E M n, W m , 6, Q M ) (4) 

where 

E M n = (yj,MVj,N - 5mn)/2 , W M = ~4>m , ©m = 6,m (5) 

are the finite strain tensor, the reference electric potential gradient, and the 
reference temperature gradient; of course, 5mn is the Kronecker delta, and 
(f) is the electric potential. Hence, by using ip the constitutive relations (4) 
of [13] are deduced for K^, A^, rj; here we rewrite them from [13]: 

dip 1 

Ku = Ui,ApoTr^ \~ JXlj E (EjEi — -EiEiOji) , 

oh, AL 1 

Al^EoJXljEj-Po-^-, (6) 

dip 

Recall that the heat-flux constitutive relation <^j 2 is restricted by 

Ql&l<0. (7) 
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Note that, in particular, (j3J) 2 includes the case in which Qm is linear in 
Ql, that is, 

Qm — —k>ml(9, Wa) ®l ■ (8) 

2.2 The initial boundary value problem for a thermo- 
electroelastic body 

To describe the corresponding boundary conditions to add to the field equa- 
tions ([I])-©, three partitions (Sa, S i2 ), i — 1, 2, 3, of the boundary surface 
S = dB can be assigned. For mechanical boundary conditions, deformation 
ft and traction ti per unit undeformed area are prescribed, respectively, 
on Su and Su; for electric boundary conditions, electric potential and 
surface-free charge A per unit undeformed area are prescribed, respectively, 
on S21 and 622; while for thermic boundary conditions, temperature 9 and 
normal heat flux Q per unit undeformed area are prescribed, respectivelyon 
531 and 532. Hence, we can write 

Ui = yi on , K Li N L = Ki on 512 ('mechanical') , (9) 

= on 521, A L iV L = -A on S 22 , ('electric') (10) 

9 = 9 on S 31 , Q L N L = Q on S 32 ('thermic'), (11) 

where N = (Nl) is the unit exterior normal on S and 

S a U S l2 = S, S il nS a = <l> (z = l, 2, 3). (12) 

We put 

A body ■= (/*, Pe,i) , (13) 

A SU rf := (ft, Ki, 0, A, 6, Q) , (14) 

A := (Abody, Asurf) = (/», PE, 1, ft, ^, A^jQJ. (15) 

•Afiod?/; Abody, and ^4. are said to be the (external) body-action, surface- action, 
and action, respectively. 

The initial conditions have the form 

j/i(X,0) = /i(X), ft(X), 0)=^(X), 
9(X, 0) = fc(X), 0(X, 0) = /(X) (X e B, t = 0) , (16) 

where 

2" = (/i, 9h h, l) 



4 



are prescribed smooth functions of domain V. The initial boundary value 
problem is then stated as: assigned Abody, to find the solution (</>, 9, yi) in 
B to the constitutive relations and field equations (Gp-(G|) which satisfies 
the boundary conditions (TP|)-( f77p and initial conditions (To}) for given A sur f 
and X. 

3 Biasing and incremental fields 

In incremental theories three configurations are distinguished: the reference, 
initial and present configuration. 

3.1 The Reference Configuration 

In the reference state the body is undeformed and free of all fields. A generic 
point at this state is denoted by X with rectangular coordinates X^. The 
mass density in the reference configuration is denoted by p a . 

3.2 The Initial Configuration 

We put 



A° := ( A body , A surf ) = (/°, p° E , 7°, Vl Kf, 4>°, A , 9°, Q° ) . (19) 



A% ody , At, ody , and A° are said to be the (external) body-action, surface- 
action, and action, respectively. 

In this state the body is deformed finitely under the action of the pre- 
scribed action A°. The position of the material point associated with X is 
given by 



In this state the electric potential, electric field and temperature field are 
denoted by 0°(X, t), W° = -<j)° a and 0°(X, t), respectively. 
The initial fields 




(17) 




(18) 




Jo = det(y° a>L ). 



y° a = y° a (X,t), 0° = 0°(X,t), 9 = 9°(X,t) 



(20) 
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satisfy the equations of nonlinear thermoelectroelasticity ffl"l)- (fT2l under the 
prescribed action A°. In studying the incremental fields the solution to the 
initial state problem is assumed known. 

3.3 The Present Configuration 

To the deformed body at the initial configuration, infinitesimal deformations, 
electric, and thermal fields are applied. The present position of the material 
point associated with X is given by ?/j(X, t), with electric potential 0(X, t) 
and temperature 6*(X, t). 

The fields yi(K, t), 0(X, t), 6*(X, t) satisfy ([I])-© under the action of 
the external action (1151) . 

3.4 Equations for the incremental fields 

Let e be a small and dimensionless number. The incremental process e{y 1 1 1 , 9 1 ) 
for (y, 0, 9) superposed to the initial process (y°, 0°, 9°) is assumed to be 
infinitesimal and, therefore, we write: 

y i = S ia (y a + ey 1 a ), = 0° + 50 1 , 9 = 9° + eO 1 , (21) 

Corresponding to ( 12TI) . the other quantities of the present state can be written 
as: 

A = A° + eA 1 , (22) 

where, due to nonlinearity, higher powers of e may arise. For the incremental 
action we have 

Ao d y--=(flp l E,l l ) (23) 

A\ urf := (yl Kl 4>\ K\~9\Q 1 ) (24) 

A 1 := (A body , A surf ) = (ft, p x m 7 1 , vl Kl 1 , A 1 , ~9\ Q 1 ) (25) 
We want to derive equations governing the incremental process 

(u:=y\ 0\ 9 l ). 
From (121jl and ( 1221) . we can further write: 

E KL = E° KL + eE 1 KL , 

Wl — W£ + eWl , (26) 

0L = e2 + ^ei, 
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where 



E° KL = (y a , K y° a ,L - Skl)/2 , E l KL = (y° a>K yl L + y° a , L yi K )/2 , 

W° L = -<j>% , Wl = -4>\ , (27) 



Substituting (|2B -( J27T) into the constitutive relations ([I])-©, with some 
very lengthy algebra, the following expression are obtained [T5] : 

K M i — 5ia(K° Ma + eK l Ma ) , A M — A M + eA M , 

rj^tf + eri 1 , Q M = Q° M + eQ l M . (28) 

where 

^Afa = GMaL-yUj,L + RLMa^L ~ Po^Mcfi , (29) 

A M = Rmn-/U 7i n — Lmn4>]n + PoPm® 1 , (30) 
i] 1 = A Mt u 7i m - Pm4>) m + a01 > ( 31 ) 

<3m = ^-MNa u a,N — Bmn4>]n + CW? 1 + F~MN@]n ■ (32) 

By putting 

K>MNa — —A MNa , K>MN = E>MN, K M = ~Cm, k MN = —F~MN, 

the latter rewrites as 

Qm = —HMNaU a ,N — k MN^]n ~ K M^ ~ ^MN^N ■ (33) 



In (1291 - (1321 . GmoL-j are the effective elastic constants, Rlmo are the ef- 
fective piezoelectric constants, Aj^ Q are the effective thermoelatic constants, 
Lmn are the effective dielectric constants, Pm are the effective pyrolectric 
constants, a is related with the specific heat. Their expressions are [15] : 



G^aL-y = Va,MPo^ (0° , E AB , W A ) y° L + p Q — (9°, E° AB , W%) 8 al + gKaL-y 



'dW^dE ML 



Rlmj = ~Po aixr aTP — {9°, E AB , W A )y° M + r KLl , 



^ = -^(e°,E AB ,w 1)yU . 
JMN = - po dw M dw N {e °> E ° AB > w ° a) + lMN ' Pm = ~dw^m {B% E ° AB > w ° a) ' 
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d ^tno TJ/0 ^ , __ 9Q M 

A MNl - — 

Bmn = Sjrf (0°, E° AB , WX) = K-mn J 



a = -W {9 °' E ° AB ' W ° a) ' AmN " = dE^ {6 °' E ° AB ' W ° a)v °> l = ~ KMN ^ 



C m = ~^f~(®°> E ABi W 2) = ~ K M , FMN = E ABi W l) = ~ K MN 

where 

9Ko.L^ = £oJo W°Wp{Xx,pXL n — Xk^X^p) + WpW°{XK )a XL,l3 — Xk, (3X1,0) 

+W° p W%{X Ktl X L>a - X K>a X L ^)/2 - WZW°X Kt pX L)f} ] , 
r KL-y = £oJo(w°X K a X L ^ — W°X Kn X L a — W°X Ka X La ) , l MN = e J X M ^ a X N) 



In (1321) we have introduced the ^-notation to allow comparison between the 
proofs written here and those written in [2] . The following symmetries hold: 

(36) 

3.5 Restriction on the incremental heat flux 

Now we show that the restriction on the heat flux (j4j) 2 , together with the 
condition 

Q° L = for 02 = 0, (37) 
implies an analogous restriction on the incremental heat flux (1321) . that is, 

Ql®l<0. (38) 

Indeed, substituting 

Q L = Q° L + eQ{, ® L = ®l + £ Q l L 

in (J7J), we obtain 

(Ql + sQ\) (Q° L + eOl) < , (39) 



which for G° L = 0, by (|37jh yields 

Note that the choice (jHJ) for the heat flux response function satisfies ( l3~7[ ). 
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3.6 Incremental field equations 

By substituting flZQ-dH into (HJ)-© and ©-(HI]), we find the governing 
equations for the incremental fields 

K lla, M + Pofa = Po '< , (40) 

*m,m = Pb, (41) 

Po (6°f^ + e 1 ?) = —Qm,m + Pol 1 - (42) 

Introducing the constitutive relations (I29 l -( l32j) into the incremental equa- 
tions of motion ( 140|) . the equation of the electric field (JHJ), and the heat 
equation (T42|) . for /* = we have 

Gmo.Isi U "i,LM + RLMa<P t LM ~ Po^Mcfi,M = PoU a , (43) 

+ PoPmOm = Pe , (44) 

= KMN^NM + k m0]m + k MN®]nM + K MNaU a ,NM + Po 7* ■ (45) 

4 Uniqueness theorem of the solution of the 
incremental differential equations 

In the present Section we assume if = and Q° L = 0, i.e. the initial 
temperature field 8° is uniform. This holds true when the initial state is 
static. We follow step by step the proof of Nowacki [2] and put in evidence 
any difference when it will appear. 

A modified version of energy balance is needed. It follows by substituting 
the virtual increments by the real increments 

5u a = — - dt = v a dt , 5u a m 
at 

in the principle of virtual work 

/ (fa-PoU a )Su a dV+ K a 5u a dS = 
Jv° K 7 Js° 

Thus the fundamental energy equation 

/ (fa - PoV a )v a dV + K a v a dS 
Jv° K 7 J s° 



u a M dt 



[ K l Ma 5u^ M dV. (46) 
Jv° 



[ Kl Ia u atM dV (47) 
Jv° 
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is obtained, where we substitute the constitutive relations f[2"9l. Hence 



J (fa - PoVo)v a dV + J K a v a dS 

= (GMaLyU-y,L + RlMo^L ~~ PoAa/q^ 1 ) U a , M dV , (48) 

thus 

j (W+K) = J jlv a dV + J o K a v a dS + J o (poAMaO 1 - Rlm^l) u a , m dV , 

(49) 

where W is the work of deformation and K, is the kinetic energy: 

W = - / G MaLj u a>M u^ L dV , JC = - f p Q v a v a dV. (50) 
2 Jv° 2 Jv° 

Now, to eliminate the term J vo p K Moi 6 l u a ^u dV , we multiply by 9 1 the 
heat-conduction equation (1451) . where if = 0, and integrate over V°; after 
simple transformations we obtain 

E 

I Po 6 1 A Ma u a> M dV =^ML [ e l cj>\N M dS + 

JV° u J S° 

+ it I o 1 N L ds+ f ^ I Po\N u ds + ^Sr7 e\ a , L N M ds (51) 

U° JS° u J S" v J S° 

+ P L Po 6 l ^ L dV + j o poO^UV - jV - (x + Xe + X<f> + Xu), 



where 

V = — , 

29° Jv 



[ poO 1 6 1 dV , (52) 
Jv° 



Xe 



X<t> 



L e ^ dv ■ x " = ^ L Ldv - (53) 

Note that this equation differs from the corresponding Eq.(25) in [2] by the 
terms x<f>, X an d Xu- Now, substituting (IBTj) into (T451) . we are lead to the 
equation 

j t (W + K + v) + ( X + Xe + Xj> + Xu) = J v J l a v a dV + J so K a v a dS + 

+ ir 1 0l ^ N - ds + K i 1 ^ ds + T I 01 ^ Nm dS + (54) 

+ ITo I Po d V dV - [ ( R LMa ^L < m - PoPm V m ) dV . 
0° Jv° jv° v ' 
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To eliminate the term 



J o (RLMa<f>l L Ucx,M ~ p o PAl0 1 (})] M )dV 

in Eq. (I54p we substitute the constitutive relations (I30p into the time-derivative 
of the equation of the electric field (14T1) with p E = . Multiplying the ob- 
tained equation by <p l and integrating over the region of the body, we obtain 

/ k M 4> x N M dV + I A M W l M dV = 0. (55) 
Js° Jv° 

Using the relations flHUj) and (|55p . after simple transformations we obtain 

/ A L w l L dv = 

Jv° 

= (R L MaU Q ,MWl + L LM W x M Wl + Po p L £(9 l Wl) - PoPlO'wI) dV = 

= - / AlN L ^dS, 
Js° 



from which 



J (R K MaU a , M Wk - Po p K e l w}c)dv = 

= - I A^dS - ±S - ±( Po P K \ y /WldV 

(56) 

where 

£ = \l km I W x M W x K dV . (57) 

L JV° 

In view of Eqs. (1541) and (1561) . we arrive at the modified energy balance 

j t (W + JC + V + £ + Po P K J v fW l K dV) + ( X + Xe + X<f> + Xu) = 

= I fiv a dV + f K a v a dS + 
Jv° Js° 

+ ^ I <Wl*m 1 1 N L dS + ^ I W L N M dS + (58) 

The energy balance (1351) makes possible the proof of the uniqueness of the 
solution. 
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We assume that two distinct solutions (u[, (fi v , 6 1 ') and (u", (ft 1 ", 6 1 ") 
satisfy Eqs. p0"]) - P"2"j) and the appropriate boundary and initial conditions. 
Their difference 

(u^u't-u'i = 1/ -0 1 ", § = V = 1 ") 

therefore satisfies the homogeneous equations (H0l - (|42l) and the homogeneous 
boundary and initial conditions. Equation (15H]) holds for (u i: 0, 9). 

In view of the homogeneity of the equations and the boundary conditions, 
the right-hand side of Eq. (|58l) vanishes. Hence 

j t (w + JC + V + S + p P K j v /W l K dv) = -( X + Xe + X<t> + Xu) < 0, 

(59) 

where the last inequality is true since by fl5*Hl) and flHHj) we have 

-(X + X0 + X4> + Xu) = j a j vo Q l M ® l M dV. (60) 

The integral in the left-hand side of Eq. (f5"9"]) vanishes at the initial instant, 
since the functions Ui, <ft, 9 satisfy the homogeneous initial conditions. On 
the other hand, by the inequality in ( |59l) the left-hand side is either negative 
or zero. 

Now we assume (i — Hi) below; note that (Hi) is the sufficient condition 
of J. Ignaczak, written in [2] on pages 176-177. 

(i) The initial deformation y° a realizes that the tensor GmuL^ is positive- 
definite, so that W > by (ED}. 

(ii) The tensor Lkn is positive-definite so that, by ( |571) . £ > 0. 

(in) Lu is a known positive-definite symmetric tensor, gj = p Q Pi is a 
vector, and c = p o a/20° > 0; consider the function 

A(6\ W L ) = (6 1 ) 2 + 26 1 g I W} + LjjW}W} 

A is nonnegative for every real pair (8 1 , Wl) , provided 

\gi\ < c\ m 

where X m is the smallest positive eigenvalue of the tensor Ljj. 
Under these three assumptions, (15*9"]) implies 

Ui, L = 0, 6 = 0, W L = 0, 

which imply the uniqueness of the solutions of the incremental thermo- 
electroelastic equations, i.e., 

u'i = u'l, e 1 ' = e 1 ", w]' = w\" . 

Moreover, from the constitutive relations we have that 

X x i — K y " A 1 ' -A 1 " n 1 ' - n 1 " 
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5 On the generalized Hamilton's principle 

We define the free energy, electric enthalpy, and potential of the heat flow 
respectively by 



Ip 1 — -GMaL- 1 U C() MU li L + RLMa4> l L U a,M — Po^ A-MaUa, M ~ Pm4> 1 M + 77 ^ 



and 



Whence 



dH 1 x dH 1 



du 



a, M 



dWl 



-a! 



L J 



dH 1 



-pov 



Q 1 --?£- 



M 



(61) 
(62) 

(63) 

(64) 
(65) 

(66) 



Lastly we define two functional 

n = J vo (h 1 + prfe 1 - flu a )dv - f (k l a u a - Ay)d£ 

and 

* = / (r-p o (r ? 1 O 1 + r ? 1 O 1 + ^ 1 1 + 7 1 ^ 1 ) N )dy + / ^QdS, (67) 
Jv° ^ ' Js° 

Eqs.dSD-dniD generalize Eqs.gl (36)-(38)]. 

The generalized Hamilton's principle has the form 



5 jf (/C - n) dt = , 8 jf * tft = 



(6* 



The virtual processes 

(<5u a , 50 1 ) 

of the body must be compatible with the conditions restricting the process 
of the body. Moreover the virtual processes must satisfy the conditions 

5u Q (x, ti) = 5u a (x, t 2 ) = 0, ^(x, tt) = ^(x, t 2 ) = 0, ^(x, t x ) = ^(x, t 2 ) = 0. 
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Hence, performing the variations in the second of Eqs. (JbSj) and observing 
5H 1 = K l Ma 5u a , M - PoV 1 ^ 1 + A l L 5$K , (69) 



that 
and 



dt 

It! 

we have 

rt 2 



t-2 



IC-Ujdt = 

o {^u a u a ~ H 1 - prfO 1 + f l a u a )dV + / (K l a u a - AY)dS], (70) 



6 ^ (K -U)dt = jT dt[J vo [- Po uJu a - KlfJu^M - A\5^ L + f*6u a )dV 



+ / (R^-AH^dS 

Hence by the identities 



(72) 



we have 



5 jf (zc-n)rft 



dt 



v° 



( - p U a Su a + K l Ma ^ M + + 



AT, Af 1 



IS" 

Thus we have 



s° 



dS 



rt2 




/ dt 




It! 





p u a + K l M M + fl)5u a dV 



v° 

+ / (K 1 a -K 1 Ma N M )5u a dS - I (A 1 + A 1 M N M )5(j) 1 dS] =0. (74) 

Since the variations 5u a and 5(f) 1 are arbitrary, Eq. (J74"]) is equivalent to 
the equations governing the incremental motion and electric field, completed 
by the appropriate boundary conditions. These equations and boundary 
conditions coincide with those written above. 

Next we perform the required variation in the second of Eqs. fl68l) by ob- 
serving that 



5T 



du. 



89 1 



l a,N w L ul f,L 



(75) 



(71) 



dV 



(73) 
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By flBTj) we have 



ft2 

5 M 



dt[ J v (st - Porfie^e 1 + e°59 l ) - p -q {e 1 5e 1 + e 1 se 1 ) - p o7 W))dv 

+ / 56 l QdS 
Js° 

I 2 dt\ ( ( - K MLa e x M bu^ L + Q\86\ - Kml9 1 mS4> X l ~ kmO^SO 1 



+ / 56 l QdS 
Js° 



(76) 



Note that 



t-2 



0, U,T=0, 1, 



since = at ti and £2- Also by using the identity 

(a L b)^ L = a LjL b + a L b iL 

we obtain 

5 M = 
Jti 



(77) 



(78 



£ 2 dt[J vo (Ql L + Potfe + f } °e 1 -jtye^v - (Q l L N L -Q)5e l ds 



t-2 



dt 



f yo (^mlJ] m ^l + k$jO)M)j + kl^ 1 ) ^](79) 



with 



KM La l M Su (x,LdV = K MLa - / 6 1 ML 5u a dV + / 9 l M N L 5u a dS ,(80) 



K E ML Q] M ^\dV = k e ml 



9 1 ML 8^ 1 dV + / ^i^ 1 ^]. (81) 

yo ' Jgo ' J 



Hence, by performing the variation f!79|) with the variations <5w a , that 
vanish, and with 56 1 arbitrary, we obtain that ([7TJ]) reduces to 



*2 



(5 / mdt = dt 

Jti Jti 



V 



s° 



Q\N L - Q)60 l dS] . (82) 
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Thus (i) the variational equation (E3j) 2 performed with 

Su a = = (83) 
is equivalent to the entropy balance 

QIl + Po^ + V^ 1 -! 1 ) = (84) 
and the boundary condition for the heat flow 

QlN L = Q, (xgS) (85) 

if and only if 

k l = . (86) 



Alternatively, by performing the variation flo^l) ^ with all the variations 
5u a , 8(f) 1 , S9 1 arbitrary, we deduce that 

(ii) the variational equation [68\) o is equivalent to the entropy balance 
[Efy an d the boundary condition for the heat flow (85\) if and only if 

k l = 0, nf 4L = 0, K MLa = . 

6 Theorem of Reciprocity of Work 

Next we extend the theorem of reciprocity of work following some steps in [2] 
on pages 179-182, where it is referred to linear thermoelectroelasticity in a 
natural configuration. Here there are some essential changes imposed by the 
presence of the initial fields. We assume that the body is homogeneous and 
moreover that the initial state is static, so that in particular 9° = 0, if = 0. 
Here we do not assume that 9° is uniform. 

The Laplace transform of functions v = z/(x, t) , 

roo 

77(x, p) = / e~ p V(x, t) dt , (87) 
J o 

will be used below. 

Consider two sets of causes A 1 , A 1 ' for incremental processes, and re- 
spective effects (u a , (ft, 9), (u' a , 0', 9'). Starting from the equations of motion 

K La,L + Pofa = PoU a , (88) 
K l La,L + PofL = PoU a , (89) 
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taking their Laplace transform, multiplying each by 0°, then multiplying the 
first by u' a and the second by u a , and making the difference of their integrals 
over the instantaneous region V, assuming that the initial conditions for the 
displacements are homogeneous, we obtain the integral equation 

J vo ¥{F a u' a - F' a u a ) dV + J o F(^7< - KKTluZ) dV = , (90) 

where F a = p a f a , F' a = p f' a - Now, by the identity (1751) and the divergence 
theorem, we have 

J vo &(KL^K ~ KKTlu^) dV = - KKu^j N L dS 

- J vo (Wa(&K),L - K\ a {e^), L ) dV , 

hence 

/ ^(l^<-KK^u^dV = [ W(l^-KEu^N L dS 

J \f° J S° 

■ j yo (&) , l - K^7aui) dV - J v &(l^(^),L-KEm,L) dV.(91) 

Hence by the latter equation and the constitutive relations ( 1251 . Eq. (l9"U]) 

becomes 

T°(F a u' a - F' a u a ) dV + T°(K\ a u' a - K l i a u a ) N L dS 

+ J ^[po^La^'u^L-tvfa^) + ^LJV 7 (^W - V^Wfy] dV 

- j vo (¥) )L (l^ a -l^u a ) dV = 0, (92) 
which is the analogue of Eq. (54) in [2] . 



Next we shall make use of the heat-conduction equation (1841) for both the 
systems of loadings, rewritten in the form 



-(^m)-M 1 = -Po{j;), (93) 

since we have 

rf = . (94) 
Hence by Eqs. flHHl) and flUTj) we obtain 



[KLNa—j^- + I^MN—^T + K L~^ + «M7V Qo 



,NM 



PP (Am 7 « 7 ,m - Pm<P] m + a 6 l ) = ~P"(~h) ■ ( 95 ) 



9 C 
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Multiplying the latter by 9° we have 



AT/ U a ^ NL E 4>,NM ®,L . 9]nM 

V [K-LNa ~Q~ Q + I^MN—^T + K L — + KMN~^~ 



Write the latter equality for both the states, multiply the first equation by 
9 l and the second by 6 1 ; we obtain 



1 ( K-LNa + K MN —^~ + K L — + K M N Qo 



-p P yeo{K Ml u^ M - P M <p] NM + ad 1 ) = -0V0o po Q_) 7 (97) 



and 



-nTn3( U a)NL ' E 4>)mN . ®)l . ® t MN 

\ K LNa 0O + K MN —^~ + K L — + K MN —^- 



-VPoO^^m^^m - Pm^+oO 1 ') = -0 1 6°Po(^) ■ (98) 

By taking the integral over V of the difference between the last two equa- 
tions, we obtain the analogue of Eq. (57) in [2] , that is, 



+ I Po0 3 (0 T '^- - ^t-) dV = 0. 
Jv° v 9° d°> 

(99) 

Finally, we make use of the equation for the electric field 

AV = , A 1 '^ = . (100) 

Multiplying both by 8°, the first by <j)', the second by <f), subtracting the 
results and integrating over the region of the body, we obtain 



yo 



& x l,l {9° cf) 1 ') - A\ L {6° <f)i))dV = 0. (.1.01: 
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By the identity (1751) we have 



6° A 1 ^ 1 '- A l ' L <p l )N L dS- 
s° ^ ' JV° 



and thus 



dV = 0, 
(102) 



S° 



9° AW - AV0 1 )N L dS - / (0°),z (A 1 ^ 1 ' - A 1 L '(j) 1 )dV 



v° 



0" 



v 



AWU " ^ V l(^),l dV = 0,(103) 



9" (Ah 

o \ 



VI' 



- A 



u ^1 



(0°) jL a 1 ^ 1 ' - Ai'^0 1 



+ J v e°[A l L Wl' - A v L Wl)dV = 0.(104) 



Now we substitute the constitutive relation 



A 1 L = -Rljv 7 w 7 , at - L LN (j) l N + p P L 9 l 
in the third integral of the last equation. We obtain 



vi' _ A x i6 1 )N L dS 



(9°) !L [Ahcf> v - A 1 'l (p 1 ) dV 



HA 1 ! 

o \ 

J vo ¥ [(r ln ^W^ - L LN ~fi~ + Po P L W)W[' - (Rln^ ~ Lln¥^ + PoPlO 17 ) W}\ dV = 0(105) 



Thus 



+ / Q° 



J 0°(& L <P' - A^^)N L dS - J vo (6°) L (Ai L 4> v - A» L p)dV 



dV = 0. (106) 



This equation is the analogue of Eq. [21 (61)]. 

Taking the expression for 

J ¥ R LNy (u^wW[ f - Rln-KnW^) dV (107) 
deduced from (jlUfip and inserting this into f[9"2"]) yields 

- J ^[po^La^u^L - r u a , L ) dV = 

J vo e~-{F a u' a - T a u a ) dv + f so e~-{K\ a u' a - Ktp a )N L dS 
+ / F(AW + A^[^)N L dS - / (^(A^-A^W 
o ¥ Po P L (¥Wt - W'Wl) dV - J vo (^U{^K - K\ a ud) dV . (108) 
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Now inserting 



m 



yields 



K LNc 



IV 



9" 



9" 



dV + . E MN J v jiW>^-Wt 1 ME) <1 v + 



Jv° v 9° 



9 1 - 



.9 1 



dV + K MN / Hqv 



V 



JVM 

9° 



9 l - 



,NM 



9" 



dV + 



+ p 



+ P J vo p o o ( - d l 'PuWh + eip M wHt) dv 

¥(F a u' a - F' a u a ) dV + T°(K l La v! a - K L ' a u a )N L dS 



+ / HAW 



"o \ 



ji/ _ A 1 



&l<P)N L dS - 



(109) 



Next in the latter equality we transform the sum of the first four integrals. 
Firstly note that by flH7|) we have 

1 = / e' pt dt = 1/p, 
Jo 

9° = 0°( x ) ¥=6°/p, 

W^T^f e ~ t,ft(x ' t)dt = 7M^' 



V 



oofgnL. -eii^)dv = — [ (9 v f -9 i f)dv. 

V 9° 9° > p Jv° v 



(110) 

Hence by these equalities and the constitutive relation for the incremental 
heat flux ( 1331) . the aforementioned sum of the four integrals equals 



-/ (9 ll Ql L -9^Q% L )dV. 
p Jv° v 



(111) 



Again by the identity 

Q>b, ML = {pb,M) t L ~ a ,Lb t M 

and the divergence theorem, the sum fillip equals 



s° 



e 1 'Q 1 L -e 1 Q lf L )N L ds 



iv 



9 1 XQ l L -9\Q^ L )dV . (112) 
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By substituting the sum of the first four integrals in Eq. (11091) by (11121) . we 
obtain 



5° 



e^Ql-e^Q^'ANLds 



0)'lQl ~ ^lQ^l) dV 



V 



+ P p M / p e°( - e lf wh + e^wHi) dv 



+p 



[ 6°(F a u' Q - F' a u a ) dV + f e°(K\ a v! a - Kl a u a )N L dS 
+ J s W(a^^ - A^[^)N L dS - J (FIl (A 1 ^ - A^0T) dV 

J vo ¥ Po p L (¥wt - e^wl) dv - J vo (¥u(k^ - WQ^) dv + 



(113) 



The latter is the final form of the theorem of reciprocity of work, containing 
all causes and effects. It generalizes Eq.[2l (62)], and reduces exactly to the 
latter in case of vanishing initial fields, that is, when the initial configuration 
is natural. 
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